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Abstract. The reaction yp — > n°n°p has been measured using the TAPS BaF2 calorimeter at the 
tagged photon facility of the Mainz Microtron accelerator in the beam energy range from threshold 
up to 820 MeV. Close to threshold, chiral perturbation theory (ChPT) predicts that this channel 
is significantly enhanced compared to double pion final states with charged pions. The strength is 
attributed dominantly to pion loops in the 2n° channel - a finding that opens new prospects for the 
test of ChPT. Our measurement is the first which is sensitive enough for a conclusive comparison 
with the ChPT calculation and is in agreement with its prediction. The data are also in agreement 
with a calculation in the unitary chiral approach. 

In the second resonance region, a recent model interpretation of new GRAAL data claimed 
a dominance of the Pi i (1440) — > gN reaction process. We present very accurate invariant mass 
distributions of 7i n Q and n°p systems, which are in contrast to the gN intermediate state and 
which show a dominance of the An intermediate state. 



INTRODUCTION 

The description of the low energy properties of the nucleon as well as the study of 
nucleon resonances remain a long-standing task of hadronic physics. The unique features 
of the 2n° channel - the strong suppression of the direct production (A Kroll-Rudermann, 
Born terms,. . . ) - open new prospects to improve the knowledge in both fields. 



Chiral Perturbation Theory 



In the low energy regime where properties of the lowest lying baryons and mesons are 
studied, an approach exploiting the approximate Goldstone boson nature of the pion has 
been developed: chiral perturbation theory (ChPT) l[lll2]. This effective field theory has 
been extended to the nucleon sector (HBChPT 1 ) . In general, it turns out, that ChPT 
is in good agreement with experiment in describing Tt—N scattering (a). From the study 
of nil production processes, complementary information to the study of the single pion 
photoproduction channels can be gained. ChPT predicts that the photoproduction 
channel is strongly enhanced due to chiral (pion) loops [6] which appear in leading (non 
vanishing) order q 3 . This is a counter-intuitive result, since in the case of single pion 
production the cross sections for charged pions are considerably larger than the ones 



ChPT is used in this paper as a synonym for HBChPT 



with neutral pions in the final state. In a calculation up to order M\, the pion loops at 
order q 3 are responsible for two thirds of the total cross section [7]. This fact makes 
this channel unique, because unlike in other channels where the loops are adding some 
contribution to the dominant tree graphs, here they are absolutely dominating. In 0], the 
following prediction for the near threshold cross section was given: 

o tot (E r ) =0.6nb(( J B r -4 /ir )/10MeV) 2 (1) 

where Ey denotes the photon beam energy and EL hr the production threshold of 308.8 

MeV. The largest resonance contribution at order M\ comes from the Pi i (1440) reso- 
nance via the N*Nnn s-wave vertex. Actually, the uncertainty of the coupling of the 
Pi i( 1440) to the s-wave nil channel was a limiting factor for the accuracy of the ChPT 
calculation 0]. Therefore, for the most extreme case of this coupling, an upper limit for 
this cross was given in addition by increasing the constant in Eq. [T]from 0.6 nb to 0.9 
nb. 

Completing the overview of theoretical calculations of the reaction yp — > n 7t°p close 
to threshold, this channel is also described in a recent version of the Gomez Tejedor-Oset 
model [8]. This model is based on a set of tree level diagrams including pions, nucleons 
and nucleonic resonances. In a recent work, particular emphasis was put on the re- 
scattering of pions in the iso-spin 1=0 channel [9]. Double pion photoproduction via the 
A Kroll-Rudermann term is not possible for the 2k° final state. In case of a n~n + Kroll- 
Rudermann term, the charged pions can re-scatter into two neutral pions generating 
dynamically a %% loop. This effect is doubling the cross section in the threshold region 
and is regarded by the authors as being a reminiscence of the explicit chiral loop effect 
described before. 

In the past, two measurements of the reaction yp — > n°n°p below 450 MeV beam 
energy have been carried out iflbT fTlll. The second experiment showed an improvement 
in statistics by almost a factor 30. Nevertheless, in the threshold region the cross section 
still suffered from large statistical uncertainties (see Fig.[l]). 



Reaction Mechanisms in the second Resonance Region 

Nucleon resonances are studied in a variety of experiments in an attempt to obtain 
information on the structure of the nucleon by comparison to quark model calculations. 
Most information has been gathered through tzN scattering and it photoproduction. A 
complementary access is the double it production where the 27T° channel turns out to be 
the most selective one. Because of the vanishing charge of the 7T°, Born terms as well as 
direct production terms (A-Kroll-Rudermann, A-pion pole) are very much suppressed. 
Previously, two measurements of the reaction yp — > n°n°p were intensively studied 
in order to extract information on nucleon resonances. The MAMI results ffllll were 
interpreted by the Valencia model and gave a strong indication for a dominance 
for the Di3(1520)^ An. In a recent paper, the GRAAL collaboration reported on a 
measurement of the 2n° channel from 650 MeV up to 1500 MeV lfl2ll . These data were 
interpreted by an extention of the Laget-Murphy model 1 1211 . Despite the bad coverage of 
the Pn(1440) resonance with the incident beam energy, the authors emphasized that the 




FIGURE 1. Left: Missing mass Mx — m p for two detected n mesons for beam energies 780-820 MeV 
MeV (gray data, dashed 27T° sim., dotted 77 sim.). Right: Total cross section for the reaction yp — ► n°n°p 
(full squares) at threshold in comparison with 1 1 1] (open circles). The prediction of the ChPT calculation 
131 is shown (solid curve) together with its upper limit (dashed curve) and the prediction of Ref . 1 9 ] (dotted 
curve). 



data could be only explained by a dominance of the Pii(1440)— ► ovV reaction process. 
We present and discuss in this paper new and very precise invariant mass distributions 
of the n°n and n°p systems. 



EXPERIMENTAL SETUP AND DATA ANALYSIS 

The reaction yp — » n n p was measured at the electron accelerator Mainz Microtron 
(MAMI) JH Il4ll u sing the Glasgow tagged photon facility lfl5l and the photon 
spectrometer TAPS [1 1 VL 1 1 8fl . The photon energy covered the range 285-820 MeV with 
an average energy resolution of 2 MeV. The TAPS detector consisted of six blocks each 
with 62 hexagonally shaped BaF2 crystals arranged in an 8 x 8 matrix and a forward wall 
with 138 BaF2 crystals arranged in a llx 14 rectangle. The six blocks were located in 
a horizontal plane around the target at angles of ±54°, ±103° and ±153° with respect 
to the beam axis. Their distance to the target was 55 cm and the distance of the forward 
wall was 60 cm. This setup covered ^40% of the full solid angle. The liquid hydrogen 
target was 10 cm long with a diameter of 3 cm. Further details of the experimental setup 
can be found in ref. 111911 . 

The fP ~^ ft°ft°P reaction channel was identified by measuring the 4-momenta of 
the two 7T° mesons, whereas the proton was not detected. The n° mesons were detected 
via their two photon decay channel and identified in a standard invariant mass analysis 
from the measured photon momenta. Events were selected, were both of the two photon 
invariant masses fulfilled simultaneously the following cut: 1 \0MeV < niyj < \50MeV. 
Furthermore, the mass Mx of a missing particle was calculated (see Fig. [T]). In case 



of the reaction yp — > tz 'tz p the missing mass Mx must be equal to the mass of the 
(undetected) proton m p . Above the 77 production threshold of 707 MeV, the r\ — > 3k° 
decay is a potential background source for the 2n° channel via events where only two 
of the three n° mesons are detected by TAPS. A Monte Carlo simulation of the 2n° and 
71 reactions using GEANT3 12011 reproduces the line shape of the measured data. A cut 
corresponding to an interval of —2.5a. . .min{+2.5o,40MeV} width of the simulated 
line shape has been applied to select the events of interest. The 77 — > 3n° background 
was estimated from these simulations to be below 2% for the highest beam energy of 
820 MeV (compare Fig. [TJ and was subtracted for the cross section determination. 
Background originating from random time coincidences between the TAPS detector 
and the tagging spectrometer was subtracted in the usual way, using events outside the 
prompt time coincidence window ifl^l . 

The cross section was deduced from the rate of the 2k° events, the number of hydro- 
gen atoms per cm 2 , the photon beam flux, the branching ratio of the n° decay into two 
photons, and the detector and analysis efficiency. The geometrical detector acceptance 
and the analysis efficiency due to cuts and thresholds were obtained using the GEANT3 
code and an event generator producing distributions of the final state particles according 
to phase space. The acceptance of the detector setup was studied by examining inde- 
pendently a grid of the four degrees of freedom for this three body reaction (azimuthal 
symmetry of the reaction was assumed). In a grid of total 1024 bins the acceptance is 
100% for the beam below 410 MeV and above greater than 95% for the energy up to 
820 MeV. The average value for the detection efficiency is 0.4%. The systematic errors 
of the efficiency determination are small, because the shape of the measured distribu- 
tions are reproduced by the simulation. The systematic errors are estimated to be 8% 
and include uncertainties of the beam flux, the atarget length and the efficiency determi- 
nation. 



RESULTS AND DISCUSSION 
Chiral Perturbation Theory 

The measured total cross section at threshold for the reaction yp — > 7r°7r°p is shown 
in Fig. [T] as a function of the incident photon beam energy. The results are in agreement 
within the rather large error bars with a previous experiment lUlll . The present data 
are compared with the prediction of the ChPT calculation [7] and is in agreement with it 
lElll . although up to 20 MeV above threshold the data are somewhat lower than the ChPT 
prediction. The ChPT prediction using the upper limit of the coupling of the Pi 1 (1440) 
to (7tn) s -wave can be excluded. In the future, this might be exploited to establish a better 
constraint on this coupling by using our result as an input. The total cross section is 
also compared to the calculation with the chiral unitary model [9] and shows a good 
agreement with this latter calculation. 
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FIGURE 2. Invariant mass of n°7i° and n°p for different bins of beam energy (full squares). The 
GRAAL data is shown by the full circles. The curves show oN phase space (dotted), An phase space 
(dashed dotted) and the model calculation H (full curve). 



Reaction Mechanisms in the second Resonance Region 

The invariant mass distributions for two beam energies are shown in Fig. They are 
compared to a An phase space and a oN phase space simulations and to the Valencia 
model calculation [8]. For the o a Breit-Wigner with a pole and a width of 800 MeV 
according to the Laget-Murphy model was assumed W22I1 . The GRAAL data around 
720 MeV beam energy agrees very well with our rn^) p data, whereas in the case of 
m^o, the agreement is worse. In the m K o p distributions, the An intermediate state dom- 
inates starting already at 600 MeV beam energy. The other phase space distributions can 
not describe the data. In the m n o n o distributions the differences between the different 
reaction processes is much less discriminative. The dominance of the An intermediate 
state in the 2n° production channel seems to be the more obvious explanation, although 
no interference effects are taken into account in this simplified comparisons. This ob- 
servation is in contradiction to the claimed oN dominance in the Laget Model IU2I1 . A 
future partial wave analysis has to clarify this discrepancy, and the presented data will 
provide strong constraints for solutions in the second resonance region. 



Summary 

In summary we have measured the total and differential cross sections for the reaction 
yp — ► n°n p. The prediction of the ChPT calculation [7] is in agreement with our mea- 
sured data I21I1 . The upper limit quoted for this prediction can be excluded. This finding 
might be exploited for a better constraint on the Pn(1440) to s-wave %% coupling. Fur- 
ther on, the cross section is also well reproduced by another calculation U, where pion 
loops are dynamically generated. 



Secondly, invariant mass distributions for n n and p% are presented. In the sec- 
ond energy region, oN phase space and hence a dominant contribution of the process 
Pi i( 1440) — > op seems to be unlikely. The differential cross sections are very well de- 
scribed with a An intermediate state. This is supported by the Valencia model, where the 
dominant contribution stems from the Di3(1520)^ An process and is in contradiction 
to the Laget model 1 1211 - In a future partial wave analysis this data could provide strong 
limitations on the resonance parameters up to the second resonance region. 
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